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^ Abstract: 
> 

We study resonant pair production of heavy particles in fully hadronic final states by 
\^ means of jet substructure techniques. We propose a new resonance tagging strategy that 

\0 smoothly interpolates between the highly boosted and fully resolved regimes, leading to 

ff^ uniform signal efficiencies and background rejection rates across a broad range of masses. 

Our method makes it possible to efficiently replace independent experimental searches, 
based on different final state topologies, with a single common analysis. As a case study, 
we apply our technique to pair production of Higgs bosons decaying into bb pairs in generic 
New Physics scenarios. We adopt as benchmark models radion and massive KK graviton 
^ production in warped extra dimensions. We find that despite the overwhelming QCD 

background, the 46 final state has enough sensitivity to provide a complementary handle 
in searches for enhanced Higgs pair production at the LHC. 
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1 Introduction 



Jets are a ubiquitous component of the LHC program, relevant for precision Standard 
Model measurements, Higgs boson characterization and Beyond the Standard Model 
searches [l|[2]. In particular, searches for New Physics in multijet events are an im- 
portant element of the LHC physics program. New resonances and contact interactions 
have been searched for by ATLAS and CMS in final states with two jets jSj-jlOj, four 



jets 11 , 12 , six jets 



13 -15 , eight jets [16] and up to ten jets for the semi-classical black 



holes searches [17j[18 . 

A challenge in searches for new phenomena in multijet final states is the prohibitively 
large QCD multijet background. A range of techniques is then required in order to identify 
particular categories of jets, making it possible to reduce this background. Among those 
that have been validated and applied to searches, one can mention 6-tagging [4,19,20 , jet 
shapes for quark/gluon and other flavour identification studies 21 22 and jet substructure 
tools 23 - 29] . Stringent constraints on a variety of new physics models have been obtained 
this way, with many more expected with the full 2011-2012 dataset. 

Prom the kinematic point of view, the most common scenario is that of a heavy reso- 
nance X produced in the s-channel which then decays back into a pair of quark or gluon 
jets. However, there is a large class of models where paired production of resonances 
dominates, processes of the form — )• X — )• 2Y — )• 4 partons, with Y being another 
massive particle. The mediator X of this production might be an exotic particle from 
a new strongly coupled sector, or a resonance from extra-dimensions, such as a massive 
graviton or a radion. The Y resonance could be either some BSM particle (sparticles 
in R-parity- violating supersymmetry, colorons [30[ , axigluons [3l]) or some SM particle 
{W, Z or Higgs) that subsequently decays into quarks and gluons. 

These generic four parton processes lead to very distinct final state signatures depend- 
ing on the interplay between the masses of the two intermediate resonances, Mx and My. 
If the mass ratio is large, Mx ^ My, the Y resonances will be produced very boosted, 
and typically the decay products of each of the two Y resonances will be collimated into 
a single /at jet. On the other hand, for Mx ~ 2My, the Y resonances will be produced 
nearly at rest, decaying into four well separated jets. Existing searches assume either the 
highly boosted or fully resolved regimes, and by doing so exclude a potentially large region 
of the New Physics parameter space. 

It is the goal of this paper to design a jet reconstruction and analysis strategy that 
can be applied simultaneously to the boosted and resolved regimes. This will be achieved 
by merging the boosted-regime strategies, based on jet substructure techniques, with a 
suitable strategy for the resolved four-jet regime, based on dijet mass pairings, together 
with a smooth interpolation between the two limits. Such a strategy has the potential to 
make the experimental searches more efficient and allow a wider range of BSM models to 
be probed within the same common analysis. 

The approach that we will present here is fully general and model-independent, assum- 
ing only that resonances are pair produced and then decay hadronically, with no constraint 
on the absolute masses: indeed, at the parton level the problem turns out to be scale inde- 
pendent, and the dynamics are completely determined by the mass ratio vm = Mxj'i'My . 
Of course additional QCD radiation and confinement on the one hand and experimental 
cuts on the other break this scale invariance, but we will see that the general qualitative 
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results are robust. 

To provide a realistic application of our technique, we will examine resonant Higgs pair 
production, recently studied as a promising probe of New Physics scenarios [32ij35j . We 
will therefore derive model independent limits on BSM resonant Higgs pair production 
in the 46 final state. We will then apply these bounds in the context of warped extra 
dimensional models, where Higgs pair production is mediated by either a spin zero (radion) 
or spin two (massive Kaluza-Klein graviton) resonance. We will show that a wide range of 
the parameter space of the radion and massive KK graviton scenarios can be covered by 
present and future LHC data, and that despite the overwhelming QCD background, the 
46 final state has enough sensitivity to provide a useful handle in searches for enhanced 
Higgs pair production at the LHC. 

The outline of this paper is as follows. We begin in Sect. [2] by introducing the general 
search strategy for pair produced resonances that can be applied simultaneously to the 
boosted and resolved regimes. In Sect. [3] we review the theoretical models for resonant 
Higgs pair production in warped extra dimensions scenarios. Then in Sect. |4] we apply the 
jet reconstruction strategy both to signal events and to the QCD multijet background, 
and explore the potential for new physics searches in the 2H — )• 46 channel. In Sect. [5] we 
conclude and outline future developments. 



2 Scale-invariant resonance tagging 

Multijet signatures have long been recognized as an important channel for Beyond the 



Standard Model searches at hadron colliders 36 . The main difficulty in these channels is 



how to tame the overwhelming QCD multijet background. Searches in multijet final states 
are commonly separated into boosted and resolved regimes. An example of the former 
arises when light partons are produced from the decay of a heavy resonance. Recently 



developed jet substructure techniques, reviewed for example in [37,38 , make it possible 
to substantially improve the discrimination power in the boosted regime. At the LHC the 
advent of jet substructure methods has made it possible to study boosted production of 
the heavy Standard Model particles, like W and Z bosons and top quarks, for -v/I, the 
centre-of-mass energy of the hard process, above 1 TeV. Searches for new physics such as 



resonant production of VV or tt 39 - 42 or searches for boosted supersymmetric particles 
and colored scalars jTT,15 have also benefited from these developments. 

In this section we introduce a general strategy for jet reconstruction designed for 
searches of pair-produced resonances in fully hadronic final states, which is simultaneously 
suitable for both the highly boosted and the fully resolved regimes and that smoothly 
interpolates between them. The generic process we are interested in is the s-channel 
production of a resonance X which then decays into a pair of resonances Y , which in turn 
each decay into a pair of light Standard Model particles, labeled z, 

pp^X ^2Y ^Az. (1) 

The ratio between the masses of the X and Y resonances will determine the degree of boost 
of the Y resonances and consequently the angular distribution of their decay products z 
that will be observed in the detector. At parton level, neglecting the mass of the final state 
particles m^, the problem is scale invariant and is characterized by a single dimensionless 
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variable, denoted by 



Mx 
2My 



(2) 



which is simply the boost factor from the Y rest frame to the X rest frame. In the highly 
boosted regime, vm — ^ oo, while in the fully resolved regime, where the intermediate 
resonances Y are produced at rest, tm = 1- Schematic diagrams for the boosted and 
resolved topologies are shown in Fig. [TJ 





y 
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y Y 
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Figure 1: Schematic diagrams for the generic process pp X ^ 2Y 4z in the boosted (left 
plot) regime, corresponding to large values of the mass ratio = Mx /2My, and in the resolved 
(right plot) regime, corresponding to small values of rj\/. 

If we assume that the heavy resonance X is produced at rest, so that the laboratory 
and center-of-mass reference frames coincide, we can parametrize the four momenta of the 
X — )• YY decay with the convention that P = {pT,x,PT,y,PL, E). We then have 



Px = {0,0,0, Mx), (3) 
Mx 

2 
Mx 



Mx 

Py^ = {Py • sin By cos (jyy, (3y • sin Oy sin (py, f3y ■ cos Oy,l) , 



Py2 = — ^ (—/3y • sin 0y cos (/)y, —/3y • sin 6'y sin (/)y, —/3y • cos 1) , 

where Yi and I2 are the two decay products of the X particle, Oy is the angle of Yi with 
respect to the beam, and (py is the azimuthal angle. The boost parameters from the 
laboratory frame to the rest frame of the Y particles are given by 



/3y = \PY\/EY = ^l-l/rl, (4) 

lY = l/y^l-/32 =rM. (5) 

As one can see, the boost of the Y particles, is independent of the absolute masses of 
X and Y, and depends only on their ratio. It is in this sense that we can consider that 
the problem at hand is scale invariant: jy does not depend on any absolute mass scale. 
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Fig. [T] suggests that depending on the value of the mass ratio vm the search strategy 
should be different. For large tm, the resonances Y will be very boosted, and thus the 
angular distances of their decay products will be small, while for low rjv/ the four final 
state particles will be well separated. Since we are mostly interested in the case in which 
the final state particles z are QCD partons, quarks or gluons, we will end up either with 
two /at jets (in the boosted regime), four well separated jets (in the resolved regime) or 
one fat and two separated jets (in the intermediate regime). Fat jets are jets for which 
the substructure pattern is unlikely to have arisen from QCD radiation. 

Given that in general we do not have information on the masses of the intermediate 
resonances, we don't know a priori in which of the two regimes we will find ourselves, and it 
would be beneficial to have a search strategy that simultaneously explores all possibilities. 
It should exhibit reasonably homogeneous efficiencies and background mistag rates for any 
value the mass ratio r^f within the physically allowed range. Below, we will present such 
a combined strategy that simultaneously explores the boosted and resolved regimes. 

In order to validate the performance of the strategy that we will propose we have 
generated events for the generic process Eq. using a toy Monte Carlo simulation. The 
heavier resonance X is assumed to be produced at rest in the laboratory frame, justified 
by the fall-off at large masses of parton luminosities |43| , and to decay into the two 
intermediate resonances Y with a homogeneous angular distribution, as if it were a spin- 
zero particle. The massless decay products of the Y resonance decay are also assumed to 
decay isotropically in the Y rest frame. In this toy simulation the possible widths of the 
intermediate resonances are neglected, as well as the masses of the final state particles m^. 

In view of the later applications to Higgs pair production, we will set My = 125 GeV 
and vary Mx in a wide range, although it should be clear that at parton level the event 
classification will depend only on the ratio rj\,f • For each of the Mx values in the range 
from 250 GeV (resolved regime, = 1) to 5 TeV (highly boosted regime, = 20), we 
have generated 50K toy MC events. 

To study the performance of the jet reconstruction strategy for a realistic collider 
environment, the parton level events from the toy Monte Carlo have been showered and 



hadronized with PythiaS 44 , version 8.170. We have done this for LHC centre of mass 
energies of 8 and 14 TeV, and we include also underlying event and multiple interactions 
with the default tune 4C of PythiaS. Initial state radiation has been modeled assuming 
that the resonance is produced in the gluon-gluon channeQ Parton and hadron-level 
events were then clustered with the anti-fcf jet algorithm [45] with radius of i? = 0.5. Such 
small radii {R = 0.5 for CMS, R = 0.4 for ATLAS) are used in most experimental multijet 
analysis. 

No additional cuts will be applied to the reconstructed jets at the parton level, so 
as to avoid introducing any explicit breaking of scale invariance. On the other hand, at 
hadron-level it becomes necessary to introduce additional kinematic cuts, which explicitly 
break scale invariance. In this section we will adopt the following set of basic kinematic 
cuts for jets in hadron-level events: 

pfJ^et > 25 GeV, W''\<5, i^T = PTjet > 100 GeV . (6) 

jets 

^This is a good approximation for the warped extra dimensions models that we will introduce as 
benchmark scenarios in the next section. 
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X > 2Y > 4Z, Toy MC, Parton Level 



X > 2Y > 4Z, Toy MC, Hadron Level, LHC 8 TeV 
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Figure 2: Left plot: the fraction of events -with a given number of reconstructed jets, as a function 
of the resonance mass ratio vm (Eq. for parton-level toy Monte Carlo events. No cuts have 
been applied to the final state particles. Right plot: the same at hadron-level, -with the basic cuts 
Eq. ^ applied. Note that at parton level the only possible topologies are t-wo-, three- and four-jet 
events, so their sum is equal to the total number of events. 



In the Ht variable the sum goes over the four leading jets of the event above the p™"^ 
cut. These loose cuts have a very limited effect on the selection efficiencies expect at the 
smallest values of rjv/jf] 

In order to identify the three different regimes, boosted, resolved, and mixed, useful 
information is provided by evaluating the fraction of events -with a given reconstructed 
jet topology. We show the relative fractions of the different jet topologies in Fig. [2] for 
both parton and hadron-level, as -well as the sum of events -with two, three and four jets. 
At parton level, by construction, events can only have between two and four jets, and we 
clearly see that four-jet events dominate at low rjv/, two-jet events in the boosted regime 
for large rjv/, with the three-jet case in between. At hadron-level both the shower and the 
kinematic cuts break scale invariance, and now we can have events with fewer than two 
jets, at very low tm, due to the basic cuts Eq. ([g]), and with more than four jets due to 
parton radiation. Note though that still between 30% and 50% of the events, depending 
on the value of of r^/ , have between two and four jets. 

Given that parton showering can significantly modify the number of jets, an event 
classification based on the number of jets is not optimal under realistic conditions. Instead, 
we use an alternative classification, based on the number of tagged jets per event, that 
is, jets that are found to have non-trivial substructure. We will proceed as follows: each 
of the two hardest anti-A;^ jets in the event is reclustered using the Cambridge/ Aachen 
algorithm ^46j with Rsj = 1.3 (where the subscript sj means sub-jet), and processed with 
the BDRS mass-drop tagger [47]. This tagger has two parameters /i and ycut- To determine 
if a jet arises from a massive object, the last step of the clustering for jet j is undone, 
giving two subjets jl and j2, with rriji > mj2; if both are significantly lighter than the 
parent jet, rriji < /x • rrij and the splitting is not too asymmetric 

—2 ^^ji,j2 > Vcut , (7) 



^ When presenting our final results in Sect. [4] we will adopt more realistic cuts, in line with those of 
typical LHC searches. 
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where ARjij2 is the angular separation between the two subjets, then j is returned as 
the tagged jet. Otherwise we replace j with jl and apply the unclustering to the new j, 
repeating the procedure until we find a subjet for which the mass-drop and asymmetry 
conditions are both satisfied. If the procedure recurses to the point where it finds a 
single-particle jet, then the jet is considered untagged. We use the values fi = 0.67 and 
ycut = 0.09, as in the original BDRS paper [47] . 

Our strategy is separated into two parts: the analysis chain, which sets the flow of 
the event classification, and the quality requirements, which determine whether a given 
topology is classified as a signal event or rejected as a background event. We discuss these 
two parts in turn. 

Analysis chain. We start by examining events with at least two jets after basic cuts. 
Summarized in the flow chart of Fig. [Sj the analysis chain depends on the number of 
mass-drop jet tags, that is, jets that have been identified by the BDRS tagger to have an 
internal structure potentially not arising from QCD radiation. 

• If the two hardest jets in the event are mass-drop tagged, we examine if these two jets 
can be identified as arising from the decay products of two boosted Y resonances. 
This is established by verifying if the two Y candidates satisfy the quality conditions 
on their mass difference and angular separation listed below, in which case the event 
is assigned to the 2-tag sample. 

• If the event has a single mass-drop tag among the two hardest jets, or if the event 
had two mass-drop tags, but was not assigned to the 2-tag sample, then we examine 
whether the event can be classified as having an underlying three-jet topology, where 
the decays of one Y resonance are collected into a single jet but not those of the 
other. Events with fewer than three jets after cuts are discarded. If there is a single 
mass-drop tag, the second Y candidate is formed by adding the four-vectors of the 
other two hardest jets in the event. If there are two mass-drop tags but the event 
has been rejected in the 2-tag category, we examine combinations whereby one of 
the tagged jets is taken to correspond to a first Y candidate, while the other tagged 
jet is assumed to be a mistag and is combined with a third jet to make up the 
second Y candidate. If the jet mass and angular quality requirements listed below 
are satisfied, the event is classified into the 1-tag sample. 

• If no mass-drop tags are found in the event, or tags have been found but the event 
has failed to be assigned in either of the above categories, we examine the possibility 
of an underlying four-jet parton kinematics. Discarding events with fewer than four 
jets passing the basic cuts, we select the jet pairing such that the combination ij 
and kl of the five hardest jets in the event leads to jet masses Mij and Mf^i that 
minimizes the difference | Mij — M^i \ , and use this pairing to reconstruct the two Y 
candidates]^ If these two candidates pass the mass and angular quality requirements 
given below, the event is classified as belonging to the 0-tag sample]^ 

^This choice, as compared to selecting only the four leading jets, improves the efficiency of the resolved 
configuration since often one gluon from large-angle initial-state radiation can have a larger pr than the 
four original 6-jets. 

^Another approach would be to determine if any pairings ij and kl of the five hardest jets in the event 
satisfy the mass and angular quality requirements. If there are multiple such pairings, then one uses the 
pairing that minimizes the difference \Mij — Mki\ to reconstruct the two Y candidates. 
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Figure 3: Flow chart summarizing the basic structure of the resonance-pair tagger algorithm. The 
quality conditions are specified in the text. 



Quality requirements. To identify the event as arising from the decay of the X 
resonance, Eq. Q, additional mass and angular quality conditions are required, which 
are essential to further suppress the QCD background. Some of these requirements are 
designed so as to apply similar conditions to both the boosted and resolved topologies. 

1. We require the masses of the two Y candidates to be the same up to a given mass 
tolerance fm, to account for experimental mass resolution, as well as mass smearing 
due to underlying event, hadronization and initial and final-state radiation: 



(myi - mY2) 



(my) 



< fm , (8) 



where (my) is the average mass of the two reconstructed Y resonances. We assume 
in this work a fixed valu^for the detector mass resolution fm of 15% [4|[6,48 . This 



requirement cannot be made too stringent otherwise a large fraction of signal events 
would be discarded. 

2. In the case in which the mass of the Y resonance is known, the masses of the two Y 
candidates must lie in a mass window around My, where the width of the window 
is determined by the mass resolution of the detector. 

My{1- fm) <mYl,mY2< My{1 + fm) ■ (9) 

Since we will be considering Higgs pair-production, we will set My = Mh = 125 
GeV in the following, though this requirement has a small impact in signal events, 
and is only relevant to suppress the QCD background. 



^In realistic analysis, the typical mass resolution depends on the mass scale and jet kinematics. The 
details depend not just on the detector, but also on the experimental jet reconstruction techniques. 
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3. The separation in rapidity of the two Y candidates must be smaller than some upper 
value, 

Ay = \yYi- yY2\ < Aymax , (10) 

motivated by the fact that for a given mass of the Yl, Y2 system, background events, 
dominated by t-channel exchange, are enhanced in the forward region, while signal 
events, dominated by s-channel exchange, tend to be more central. We will take 
Aymax = 1-3 in the following, a value optimized from the high mass dijet searches 
at the LHC lillel. 



4. Likewise, the separation in rapidity between the two jets of a y candidate in the 
resolved case, yn,! and yyj,2) with i = 1,2, must also be smaller than some upper 
value, possibly different from before. 

Ay = \yYi,i - yviM < ^ylZ^ , (n) 

since for these kind of topologies, signal events will be produced closer in rapidity 
that QCD multijet production. We will take AyJ^|^ = 1.5 in the following, and we 
discuss below the rationale for this choice. 



5. To prevent excessively asymmetric configurations, whenever we have two resolved 
jets that correspond to a given y-candidate, one with p!^^ and the other with < 
p^^ (in either the 1-tag sample or the 0-tag sample), we require 

P? > ycut • P? ■ (12) 

This cut plays a similar role as the asymmetry requirement in the BDRS mass-drop 
tagger, Eq. Q, but now in the case of resolved jets, and it helps reject events where 
a soft jet arises from final-state radiation (FSR)|^ 

6. With a similar motivation, for each two resolved jets in a y-candidate with mass 
niYi, with i = 1,2 we impose the mass-drop condition on the masses of these two 
resolved jets, myj^i and mYi,2-, as follows 

max {niYi,!, mYi,2) < IJ- ■ mYi , (13) 

where fj, is the same parameter as in the BDRS mass-drop tagger. Together with the 
asymmetry condition above, applying the mass drop requirement also in the resolved 
jets ensure that the same conditions hold for the three different possible topologies, 
from the highly boosted to the fully resolved regimes. Note, however, that for our 
default choice of = 0.67, the mass-drop cut has only very limited impact on the 
final reconstruction efficiency. 

The values of the parameters used in our implementation of the jet reconstruction strategy 
are summarized in Table [TJ 

In addition, in order to improve on resolution, jet masses are filtered |47| as fol- 
lows: the constituents of each tagged jet are reclustered with a smaller radius Rmt = 



® To see the equivalence with the cut of Eq. |7|, note that for reasonably small AR, the mass of the Y 
candidate is ~ p^'py ' A7?i2. Making use of the fact that p^' < p!^\ Eq. til reduces to p^' /Pt' > Vent- 
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Jet Reconstruction 



R 


Rsi 


Ri 






Vcut 


0.5 


1.3 


0.3 


3 


0.67 


0.09 



Basic cuts 



„min 

Pt 


1 Umax 1 


trmm 

tlrp 


25 GeV 


5.0 


100 GeV 



Quality requirements 



My 


^ymax 


^ res 


fm 


125 GeV 


1.3 


1.5 


0.15 



Table 1: Upper table: parameters that define the jet reconstruction strategy, including mass-drop 
and filtering. Middle table: basic jet and event selection cuts. Lower table: parameters of the 
quality requirements imposed on the tagged resonances. See text for description. 



min (Ai?sj,sj/2, i?/), with Rf = 0.3 and Ai?sj,sj the angular distance between the two sub- 
jets after mass-drop in the boosted case. Then only the three hardest subjets, ngu = 3, 
are retained to account for at least one QCD emission. The filtering procedure improves 



mass resolution of the reconstructed resonances [47 49 and makes the procedure more 
resilient to soft radiation from the underlying event and pile- up [sOj 



This jet reconstruction strategy has been implemented in a code based on Fast Jet3 52 
and we have processed the parton and hadron-level toy Monte Carlo events through it. 
We show in Fig. [4] the efficiency of the resonance pair tagging algorithm as a function of 
resonances mass ratio vm for the parton and hadron-level toy Monte Carlo events. We 
show both the total efficiency and the breakup of the efficiencies corresponding to the 
2-tag, 1-tag and O-tag samples. The impact of the moderately loose selection cuts Eq. ([6]) 
on the parton-level efficiencies is negligible, and thus the differences between parton and 
hadron-level arise from initial and final-state radiation. 

At parton level, at low rj\,/, the O-tag sample dominates as expected from the resolved 
regime, while for large tm, the boosted regime, it is indeed the 2-tag sample that domi- 
nates. The 1-tag sample is important at intermediate boosts. The combined efficiency is 
found to be rather flat in all the mass range, between 30% and 40% for all mass values, 
showing that we are able to obtain a reasonable tagging efficiency irrespectively of the 
degree of boost of the X resonance decay products. At hadron-level efficiencies are some- 
what lower due to additional parton radiation and underlying event at low masses, but 
still we obtain a reasonable tagging efficiency of between 20% and 30% in all the relevant 
range, approximately constant for all topologies, except close to tm = 1- The 1-tag sample 
and the low r^/ O-tag sample are the ones most affected by the transition from partons to 
hadrons. 

Let us mention that the production threshold region close to r a/ ~ 1 is challenging 
from the jet reconstruction point of view. First of all, there will be a substantial degree 

^In principle for tlie resolved configuration, one could consider supplementing the analysis chain with 
the inclusion of a large-angle radiation recovery procedure, to improve on mass resolution, as advocated 

Such large-angle radiation recovery procedure leads to the so called wide-jets 
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X > 2Y > 4Z, Toy MC, Parton Level 



X > 2Y > 4Z, Toy MC, Hadron Level, LHC 8 TeV 
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Figure 4: Left plot: The efficiency of the resonance pair tagging algorithm as a function of 
resonances mass ratio tm Eq. ^ for parton-level toy Monte Carlo events. We sho-w both the total 
efficiency and the break-up for the 2-tag, 1-tag and 0-tag samples. No cuts have been applied 
to the final state particles. Right plot: same for hadron-level events, -which include the basic jet 
selection cuts. 



of overlap bet-ween the decays products of the t-wo Higgs bosons, since the two are at 
rest, which leads to -wrong mass pairings. Second, it is quite frequent that large-angle 
initial-state radiation (ISR) appears as additional jets, again confusing the pairing of the 
original jets. 

To help understand the efficiencies that we find, let us recall that the asymmetry cut 
Ucut in the BDRS mass-drop tagger leads to an upper bound on the efficiencies of signal 
events in the boosted regime of approximately ~ (1 — 2?/cut/(l + Ucut)) (this result is exact 
for the two-prong decay of a highly boosted scalar). For the two tag sample at large rjv/, 
we therefore expect that at parton level the tagging efficiency is given by 



Jim _^ . ^j/cut V exp(Aj/tnax) " 1 . , 

for our choice of parameters, where the last factor accounts for the contribution to the total 
efficiency from the cut in Ay between the two ^-candidate jetsj^ This is exactly what 
is obtained in Fig. |4j At hadron-level the efficiency in the boosted regime is somewhat 
smaller due to the contamination from initial-state radiation and the underlying event. 

When we have a resolved Y resonance candidate, the separation in rapidity between 
the two jets must be smaller than some upper value, Ay^l^. To determine the value of 
this cut, we note that in the small R limit, if such a cut is the only one applied to the 
final state, the efficiency of the fully resolved case for rjv/ ~ 1 is given by 



€0-tag {rM ~ 1) 



exp(Aygl^) - 1 
exp(Ay-|J + 1 



(15) 



so demanding that the efficiencies at low vm match the asymptotic large value, Eq. ( 14 ) 



The distribution of Y resonances is flat in cos^*, the decay angle of Y in the rest frame of the X 
resonance, thus the cut in Aj/max directly determines the maximum value of cos 6'^^^^ that wifl lead to 
resonance tagging. 
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Figure 5: The efficiency of the resonance pair tagging algorithm as a function of the resonance 
mass ratio vm, Eq. ([2]), for toy Monte Carlo events, comparing parton-level and hadron-level results 
for LHC 8 TeV. We show the total efficiency (upper left plot) and then the break-up for the 2-tag, 
1-tag and 0-tag samples. The parton level results including the basic selection cuts Eq. ^ are 
very close to the inclusive parton level case and thus not shown here. 



we obtain 

Ay-L = In ^ 1.5 (i6) 

for the default value of t/cut used in the mass-drop tagger algorithm. With this choice, 
we can achieve at low r^f the same efficiency as at large tm, at least in the parton-level 
case without the basic kinematical cuts. Note that this cut ensures not just uniform signal 
efficiency, but it also is useful for background rejection, especially in scenarios where My/2 
is substantially larger than the jet p™™ cut. 

A final interesting comparison is that of the efficiencies between parton and hadron- 
level, to gauge the robustness of our event classification based on giving priority to the 
mass-drop tags. This is useful in order to understand the impact of parton showering 
and underlying event, as well as of the basic kinematic cuts in the tagging of the heavy 
resonances. Let us recall that the only difference in the analysis chain between parton 
and hadron level events are the basic cuts in Eq. ([g]), and the fact that in the 0-tag case 
we study the mass pairings of the five leading jets. Results are shown in Fig. [5] We show 
hadron-level results only for 8 TeV, those at 14 TeV are very similar. The efficiency for 
the two tag sample is very similar at parton and hadron-level, for all the values for rjv/. 
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Figure 6: The total tagging efRciency for parton-level events (left plot) and for hadron-level events 
(right plot) for different values of the jet radii i? as a function of the mass ratio tm- The default 
radius used in this paper is i? = 0.5. 



The efficiencies for and 1-tags sample fiave a similar sliape but a smaller magnitude, 
and the shape is somewhat shifted down to lower vm values. This small shift between 
the parton-level and hadron-level efficiencies is perhaps attributable to transverse boosts 
induced by initial-state radiation. 

Since the jet reconstruction strategy that we advocate is approximately scale invariant, 
one also expects the results to be reasonably independent of the jet radius R used in the 
jet clustering: while the relative fraction of 2-tags, 1-tag and 0-tag events will of course 
vary with R, their sum should be stable. Indeed, we show in Fig. [6] the total efficiency in 
parton and hadron-level events for three different radii, R = 0.3, 0.5 and 0.8. At parton 
level, results in the boosted regime are strictly i?-independent, as shown in Fig.|6j Except 
at very low masses, parton-level results are i?-independent in all the mass range. For low 
rjv/, the degradation at parton level with increasing R arises in part because the likelihood 
that the decay products from different Y resonances end up in a single jet is higher for 
larger R. Also at hadron-level the total tagging efRciency is reasonably independent of R. 

To summarize, in this section we have presented our general strategy for a resonance 
reconstruction analysis that can be applied simultaneously to the boosted and the resolved 
regimes, with a smooth transition between the two limits. It is clear however that some of 
the details of the strategy can be modified without affecting the general philosophy. One 
could study different ways of dealing with the four-jet events rather than selecting the 
pairings which minimize the relative dijet masses, like cuts in the angular distributions. It 
is also possible to extend the number of jets considered to build the resolved F-candidate 
in the 1-tag case up to the fourth or the fifth jet, in analogy with the procedure used for 
the 0-tag case. 

These modifications could lead to an overall improvement of the tagging efficiency, but 
the basic strategy would be left unaffected. Finally, other substructure taggers could be 
used to classify events, such as N-subjettiness 53 or pruning 54 among many others 



(see [55] for a recent systematic comparison). Note however that those taggers with an 
asymmetry cut, like mass-drop and pruning, are special, because that cut can be linked 
with 3- and 4-jet analysis parameters, as done in the present analysis. In this respect, 
N-subjettiness is quite different, because it is cutting on the radiation pattern in the jet. 
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3 Resonant Higgs pair production in warped extra dimen- 
sional models 



Now we discuss the benchmark models that we will consider for resonant Higgs boson pair 
production. These models are based on the warped extra dimensions scenario |56| , where 
Higgs pair production is mediated by either a Kaluza-Klein (KK) graviton or by a radion. 
We will assume that the Higgs is the Standard Model boson f57] , and consider its dominant 
decay into two pairs of bb quarks. Higgs pair production in the Standard Model has a small 
cross section |58] (approximately 18 fb at 14 TeV), but larger rates can be expected in 



New Physics models [34,35,59] like supersymmetry, composite models, and warped extra 
dimensions. With this motivation, in this section we review the theoretical expectations 
for resonant Higgs pair production in the context of warped extra dimensional models, 
keeping in mind that the strategy proposed in this paper is equally valid for any other 
Higgs pair production scenario. 

Due to Bose symmetry, only resonances of spin zero and spin two can decay on-shell into 
a pair of Higgs bosons. Both types are present in models with warped extra dimensions. 
They are referred to as radion and KK-graviton, denoted by (j) and G respectively. These 
models can naturally explain the large hierarchy between the Planck and electroweak scales 
by introducing a nontrivial geometry in the extra dimension. The background metric for 
the case of a single extra dimension is given by 

ds'^ = e-^''yr]^^dx''dx'' - dy\ (17) 

where y refers to the coordinate in the 5th dimension and k is related to its curvature. The 
so-called ultraviolet (UV) and infrared (IR) branes are introduced at y = and y = L, 
respectively. Depending on the scenario, SM fields can be localized in the IR brane or be 
allowed to explore the 5th dimension as well. At each position in the extra-dimension y^, 
the local cutoff is given by [60] 

A{y.) = e-''y*A{y = 0) . (18) 

If one assumes that the fundamental scale of the theory A(y = 0) is the reduced Planck 
mass Mp ~ 2.4 x 10^® GeV, this scale would be locally shifted to lower values as one 
moves into the extra-dimension. Physics at the IR brane would have a cutoff at the TeV 
scale, if one requires a mild tuning of kL = 0(50). This is the solution to the hierarchy 
problem in warped extra-dimensions, where Planck-scale physics appears as TeV physics 
via the warping in the local cutoff. 

As in any model with compact extra dimensions, one expects the existence of a tower 
of massive resonances for each particle propagating in the extra-dimensions, called the 
Kaluza-Klein (KK) resonances. In particular, the graviton and the radion (and their 
associated KK towers) are described by the tensor and scalar quantum fluctuations of the 



metric, introduced as an expansion in Eq. (17) 



The fluctuation of the size of the extra dimension y is described by the 4D scalar radion 
field, denoted here by (p: 

F{x,y)(xe^''y^{x) , (20) 
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Profiles 




J(y-L) 



Figure 7: Relative localization of the graviton (blue-dashed) and radion (red-solid) in the extra- 
dimension. The profiles shown are the y-dimcnsion component of the wave functions Eqs. 20 
and[2il 



where (l){x) is the 4D wave function and e^^^ the localization profile in the fifth dimension. 
The fluctuations of 4D space-time is described by the graviton field G^^(x, y). The massless 
zero mode of this field corresponds to the usual graviton. The first massive excitation, the 
lightest KK-graviton which will focus on, is 



G^^K^,y)o.^^yj,{e^^y'^)G^^{x), (21) 



where J2 is a Bessel function and ma the graviton mass. The mass of the KK-graviton is 
related to k/Mp and to the ultraviolet mass scale of the theory Aq by 

mo = ^^^rAo, (22) 

where xi = 3.83 is the first zero of the Bessel function Ji. The three parameters, k, the 
UV mass scale Aq and the reduced Planck mass Mp are related by Ac = e~^^Mp. The 
mass scale is expected to take a value in the few TeV range. 

We neglect the effect of localized kinetic terms, which would change the value of xi [61| . 
We show the radion/graviton localization profiles in Fig. [7j both fields are localized to- 
wards the infrared brane, but the graviton localization is stronger. 

The radion and KK-graviton couplings to the SM particles are fixed by the action 

S = -\j d^xdy^6gMNT^'^ (23) 

where M, N are 5 D indexes and T^''^^ is the 5D energy-momentum tensor involving all 
fields. After dimensional reduction, the effective coupling between the radion and KK- 
graviton lightest mode and the SM is given by: 

£ = - ^(pTH' (24) 

Ag ^ A^ ^ 

where the T^^ are the four- dimensional energy-momentum tensors of the Standard Model 
species i = b, f,V, H, and V denotes a generic gauge boson. Here we are neglecting 
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corrections depending on the fermion localization parameters cl and cn, that are small 
when the fermions are localized in IR brane such as the top quark and they are not large 
for the b quark [62^. It is important to notice that the radion couples to the trace of the 
energy- momentum tensor, which vanishes at the classical level for massless gauge fields. 
Note also that the radion scale is related to the KK-graviton scale by = ^/GAg 



63 



The coefficients Cj and di are proportional to the wavefunction overlap of the gravi- 
ton/radion and the SM fields. For example, the Higgs is IR-localized, as its vacuum 
expectation value and mass are IR effects. The fact that the KK-graviton wave-function, 
Eq. ( |21[ ), is more peaked towards the IR brane than the radion, Eq. (20), as shown in 



Fig. [7j translates into a stronger coupling of the graviton to the Higgs, beyond what is 
expected from the trivial rescaling of and A^. 

In the original Randall-Sundrum model (RSI) all SM fields are localized in the IR 
brane, so all couplings are Cj ~ C(l)- More realistic models, consistent with experimental 
constraints, must have SM fields in the bulk, leading to different values of the couplings 



Cj. A well motivated configuration, which we will refer to as bulk RS 64 , predicts that 
the SM fields communicating to the EWSB sector are peaked towards the IR brane. This 
is the case of the Higgs and longitudinal W and Z bosons, and possibly the top quarks. 
Light fermions would be localized near the UV brane, whereas massless gauge bosons are 
de-localized. The graviton and radion would then couple preferentially to IR localized 
fields, namely h, Wl, Zl and possibly t as well. The coupling to 7 and g is suppressed by 
a volume factor ~ 1/kL, and the coupling to light fermions (including the quarks in the 
proton) would be extremely suppressed. In summary, the two scenarios we are going to 
consider, RSI and Bulk RS are defined by the following hierarchy of couplings of the KK 
graviton to SM particles: 



RSI scenario: 


ch = all the other q ~ 0(1) 


(25) 


Bulk RS scenario: 


CH ^ cz,w,t ^ 0{1) ~ (kL) c^^g > Cu,d,t.. 


(26) 



As an example, the stress-energy tensor for the Higgs field is given by 

r/^ = df^HduH - ^g^, {daHd'^H - mlH^) (27) 
which results in the following couplings, 

C D -^G^^} df^Hd'^H + ^cb {daHd"H - 2mjjH'^) (28) 

We do not consider a nonminimal Higgs coupling to gravity, which would require the 
use of an "improved" energy-momentum tensor ^65j and would lead to a Higgs-radion 
mixing 



66,67 , since mixing is not relevant in our case where we require > 2Mh- 
The Feynman diagrams relevant for Higgs pair production mediated by a radion (j) and a 
KK-graviton G^i^ are schematically shown in Fig. [s 
In the remainder of this section we will discuss t : 
graviton in the RSI and Bulk-RS scenarios. It is beyond the scope of this paper to review 
the experimental constraints on the parameter space of these models. A discussion of the 



le production rates of the radion and 
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Figure 8: Feynman diagrams for Higgs pair production in warped extra dimensions models medi- 
ated by a radion (f) (left plot) and a KK-graviton G^^ (right plot). The Higgs bosons then decay 
into a pair of b quarks. 



implications of recent measurements for limits on extra dimension models can be found 
in 23,68-70]. Note than when experimental limits arise from decays to photons, leptons 
and four-fermion operators involving light fermions, they can be interpreted only in the 
context of RSI but not in the bulk RS scenario. Bounds on the radion mass as a function 
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of the cutoff scale have been compiled in Refs. 
3.1 Production rates at the LHC 

The production rates of the radion and the graviton at hadron colliders will depend on 
the respective couplings to the light quarks and gluons in the incoming protons. In RSI, 
the graviton couples to light quarks and gluons with the same coefficient 0(1), whereas 
in the Bulk-RS scenario the couplings to light quarks are very suppressed. The coupling 
of the graviton to gluons in the bulk RS model is given by 

2(1 - Jo(xi)) , , 

whereas in RSI Cg = 1. 

The coupling of the radion cj) to gluons (and to photons) vanishes at tree level due 
to classical scale invariance. At 1-loop level it arises due to the trace anomaly, which is 
related to the beta function, and the top quark triangle diagram. We denote by Kg the 
coupling of the radion to gluons defined by 



where nt is given by 



62 



87r 4kL' ^-^^^ 
where we have neglected the top loop contributions. The coefficient of the QCD /3 function 
is 63 = 8. The RSI case corresponds to neglecting the volume suppressed term. Note that 
as compared to the graviton, for the radion production cross section there is less model 
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Figure 9: (Left plot) Production cross sections as a function of the graviton and radion masses 
(tog^a), where the trivial dependence on the scale A — Acfl and the coupling to gluons {c = Cg, 
1) is factored out. Solid (dashed ) lines correspond to 8 (14) TeV. (Right plot) Graviton cross 
sections (at 8 TeV ), for specific choices of Cg=l, 1/kL corresponding to RSI and bulk RS. 



flexibility in that the coupling to the gluons is fixed independenttly of the localization of 
the SM fields. 

In the following we will assume that gluon fusion is the dominant process for both 
radion and KK-graviton production at the LHC. This is certainly true for the radion, and 
also for the graviton in the bulk RS scenario, where the couplings to the light quarks in 
the proton are very suppressed. In the narrow width approximation, the production cross 
section via gluon fusion of a particle X with mass M is given by 

app^x{M, s) = Lgg{f)a{gg ^ X){fs)df, (32) 

where here r = M'^/s and Lgg is the gluon luminosity function. 

We computed the production cross section of the processes pp ^ G, (p at leading order 
(LO) using Madgraph5 |73]. The results are shown in Fig. [9j We plot the production cross 
section, where we have factored out the trivial dependence on the coupling to gluons and 
the scale of dimension- five operators, Aq and we also show the cross section for specific 
values of those parameters. 

Notice that the KK-graviton production cross section is larger than the corresponding 
radion cross section due to the fact that the radion coupling to gluons is loop induced, 
whereas the KK-graviton has tree-level couplings to gluons. Also, the KK-graviton has 
five degrees of freedom, compared to the single degree of freedom of the radion. 



3.2 Graviton and radion decays 

In RSI, with all the SM fields localized on the IR brane, a heavy graviton would decay 
democratically to all degrees of freedom. In the bulk RS, the Higgs and fields associated 
with EWSB are still IR localized, and using the equivalence theorem, one can show that |74| 

r(G ^ HH) = T{G ^ ZlZl) = r(G ^ W+WE)/2 = ^ ^ (33) 
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In bulk RS, the width to gluons and photons is suppressed by the effective volume kL, 



r(G ^ gg) = 8r(G ^ 77) ~ 8 



r(G ^ HH) 



10"^r(G ^ HH) 



(34) 



The graviton would also couple to fermions localized near the IR brane. In many mod- 
els, third generation quarks are pushed towards the IR brane via a localization parameter 
which is a ratio of a 5D mass term, Mf, and the curvature, v = Mj/k. The effect of 
v is as follows: for u = 1/2, the conformal value, the fermion zero mode is delocalized in 
the extra dimension, as the profile is flat and does not prefer a particular location inside 
the extra-dimension. For z/ > 1/2, the fermion zero-mode will be localized towards the IR 
brane, whereas for < 1/2, the localization is near the UV brane. 

The width to tops is given by 



r(G ^ tt) = 

for rriQ ^ 2mt. We have defined 

3 l + 2iyt 



1 f r,. ^2 



(35) 



dyy 



2 1 _ e-fcL{i+2;.t) J2(3.83) 
The branching ratio of the graviton to the Higgs depends on the top localization as 

1 1 



(36) 



BR{G ^ HH){ut) 



4 1 + fiiyt) 



(37) 



In Fig. 10 we show that the maximal branching ratio to a Higgs boson pair is 25% 
(when the Higgs is 1/4 of the IR degrees of freedom) , and quickly decreases as one increases 
the branching ratio to top quarks (increasing the value of ft, and therefore the localization 
towards the IR brane). 
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Figure 10: Branching ratios of the graviton to the Higgs (black), tops (blue) and vector bosons 
(red-dashed), as a function of the top localization parameter Vf 
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The dominant decay modes of the radion are into pairs of massive gauge bosons, Higgs 
bosons and top quarks. Since the couphngs are determined by the masses of the final state 
particles, and these masses arise from the TeV localized Higgs boson, the RSI and bulk 
RS couplings are the same at leading order. The corresponding widths (for large m,^) are: 



1 



HH) = r(</. ^ ZZ) = r(</. ^ WW)/2 = ^ ^ (38) 



and 



r(0^tt) = ^f^V^. (39) 

Hence, for large radion masses the branching fraction to a pair of Higgs bosons is approx- 
imately 25%, independent of A^, since the contribution from decays to top quarks can be 
neglected. We note that for the smaller values that are relevant for phenomenology 
the decay into top quarks should in principle be taken into account, but in this work for 
simplicity we will assume that BR {cp — )• HH) = 25% independent of the radion mass. 

3.3 Composite duals and model dependence 

So far, we have described higgs pair production via gluon fusion into a radion or KK- 
graviton in warped extra-dimensions. In this context, the graviton cross sections are 
larger than the radion by at least an order of magnitude, and there is little room for 
changing this hierarchy. 

To test how robust this prediction is, we would like to approach this model building 
in extra-dimensions from the point of view of holography. In this approach, models in 
warped extra dimensions are an analogue computer for strong interactions. This duality 
between 4D strongly-coupled theories and 5D weakly-coupled theories with gravity was 
inspired by the AdS/CFT correspondence, but took hold on a more qualitative basis |75[ 



and has been used to build models of QCD 76 , technicolor 77 , composite Higgs (78], 
and even condensed matter systems (79|, with some success. In this context, the KK 
resonances, a consequence of compactification, are the holograms of massive resonances 
due to confinement. 

The KK graviton is therefore the dual of a spin-two bound state in a strongly coupled 
theory, very much like the /2 of QCD fSO). One could then wonder how different the 
coupling structure of the /2-like resonance would be with respect to the KK-graviton. As 



was shown in Ref. 81 , the couplings we propose in Eq. (24) saturate the possibilities, 
once Lorentz, gauge and CP invariance are assumed. No other structures are allowed up 
to dimension-six operators. Hence, our KK graviton analysis can be directly generalized to 
strongly coupled sectors with spin-two resonances. Moreover, if the new strongly coupled 
sector participates in the electroweak symmetry breaking mechanism, a sizeable coupling 
to Higgses would be expected. 

The dual of the radion would be the dilaton, the Goldstone boson of scale invariancej^ 
The dilaton couplings at tree level would be perfectly mimicked by the radion couplings, 

^It is unclear whether in QCD one would have such a creature, but some proposals are the /(975) 
resonance [82] of the a particle 83 . 
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as the dilaton will couple to the trace of the stress tensor. This can be shown by writing 
down an effective theory where the dilaton is spurion of the scale symmetry |84| . Within 
this analysis, the loop contributions to the dilaton to massless gauge bosons will follow 
the same structure as the anomalies written in Eqs. (31). Therefore, our analysis of the 
radion couplings is also applicable to a dilaton in a composite sector. 

In summary, the structure of couplings we describe for the KK-gravitons and radion 
would be the same for the bound state duals. The main difference between the analysis in 
extra-dimensions and composite theories is the strength of the coupling. For example, one 
could imagine a composite theory where the spin-two resonance is made up of colorless 
techni-quarks, hence there would be no tree-level coupling to gluons, whereas the dilaton 
couplings would be determined by the scale symmetry. In this case, one would expect a 
larger production of dilaton than spin-two resonances. 



4 New Physics searches in the HH — > 46 final state 

In this section we apply the general resonance tagging strategy presented in Sect. [2] to 
a particular scenario, namely the resonant Higgs boson pair production pair which then 
decays into four 6-quarks. The results presented here are model independent and can 
be applied to any generic BSM model with enhanced Higgs pair production [33 
though we will provide an explicit interpretation of exclusion limits in terms of the radion 
and graviton couplings in the warped extra dimension models of Sect. [3} 

First of all, we discuss the Monte Carlo event generation for the signal and background 
events with MadGraph and Pythia and evaluate the tagging efficiency as a function of rjvf , 
to compare with the approximate kinematics of the toy MC used in Sect. [2] We recall 
that the main differences between the toy MC and MadGraph are that the latter includes 
the rapidity distribution for the X resonance and the correct angular distributions of a 
spin-two particle in the case of the graviton. Then we present the different assumptions 
that underlie our implementation of 6-tagging. A discussion of the background rejection 
capabilities of the tagging algorithm follows, where we show that the combination of the 
resonance tagging and 6-tagging reduces the QCD multijet background by several orders 
of magnitude. In the last part of the section we present the implications in terms of model 
independent searches in the HH — t- 46 final state, and interpret these results in terms of 
exclusion ranges in the parameter space of warped extra dimension models. 

4.1 Monte Carlo signal event generation 

Our benchmark model is s-channel Higgs boson pair production mediated by a radion or a 
massive Kaluza-Klein graviton resonance in scenarios with warped extra dimensions. We 
have implemented these scenarios in the MadgraphS Monte Carlo program [73]. While the 
main motivation to study both radion and graviton simultaneously is to cover a wider range 
of the model parameter space, a useful by-product is to validate the jet finding strategy 
for two different angular distributions of the decay products. Indeed, from the kinematic 
point of view the radion and graviton cases are identical (for equal masses) except for the 
different angular decay distributions of spin-zero and spin-two particles. Note that in the 
radion case, since the radion is an scalar, the kinematics and angular distribution will be 
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very close to those of the toy Monte Carlo of Sec. [2] used to validate the resonance tagging 
algorithm, with the only difference arising from the rapidity distributions of the radion. 

We have followed Ref. |85| to model the radion couplings to the Higgs boson and to 
gluons using the FeynRules framework [86] . The implementation of the model has been 
based on the default MadGraphS model with effective theory coupling of the Higgs to gauge 
bosons. In addition to the SM parameters, in the radion model we have four additional 
parameters: the radion mass, M^, the ultraviolet mass scale of the theory, A^, the radion- 
Higgs mixing parameter, ^, and the compactification scale, kL. These parameters take 
the value = 3 TeV and kL = 35, supplemented by the no mixing condition that reads 
^ = 0. The absence of mixing is justified by the fact that the radion masses considered 
will be always much larger than the Higgs mass. Any modification of A,^ translates into a 
trivial rescaling of the total rates. 

To simulate the graviton production we have used the standard Randall-Sundrum 
model as implemented in MadGraphS. Here the relevant additional parameters are only 
the graviton mass Mq and the ultraviolet mass scale of the theory, Aq, chosen to be 
Ag = 3 TeV. As pointed out on the Sect. [3| the mass scales A^j, and Aq of the radion 
and the graviton are theoretically related. However, from the practical point of view we 
select independently the parameters of the two models in order not to impose additional 
constraints on the search ranges. For both radion and graviton event generation the narrow 
width approximation has been assumed. 

We have generated events for radion and graviton production for a range of masses 
between 250 GeV and 3 TeV. Higher masses lead to too small cross sections to be of 
any phenomenological interest. As in the case of the toy Monte Carlo events, MadgraphS 
parton level events were showered and hadronized using PythiaS with the same settings 
for underlying event and multiple interactions. 

We have already discussed in Sect. [2] the tagging efficiency of the algorithm for the toy 
Monte Carlo kinematics, both at parton and at hadron level. However, the basic selection 
cuts did not match those of a realistic experimental analysis. We will use the following 
selection cuts instead in the rest of this paper: 

pf"" > 50 GeV, |?7j^*| < 2.5 , Ht = ^t^t > 300 GeV . (40) 

jets 

These cuts are inspired by typical trigger and angular acceptances of the LHC experi- 
ments 
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We show in Fig. 11 the comparison between the hadron- level tagging efficiencies at 
LHC 8 TeV between the toy Monte Carlo events and the MadGraph radion and graviton 
events, as a function of the mass ratio tm- As we can see, the toy MC results agree well 
with the radion events, which is a non trivial cross-check that event generation is under 
control. Also, the efficiencies for the radion and graviton are very similar, showing that 
the spin-zero vs. spin-two angular distributions do not lead to any large differences at 
the level of the reconstruction]^ Note that we have generated fewer mass points with 
MadGraphS than with the toy MC, hence the somewhat less smooth distributions in the 
former case. 



^''At first sight tiiis may appear to be surprising, given that the radion and graviton angular decay 
distributions are substantially difTerent. However our Aj/max = 1-3 cut is sufficiently large that the integral 
of the Ay distribution up to Ai/max is not too dissimilar in the two cases. 
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One significant difference between Fig. 11 and the results of Sec. [2| is the much lower 
efficiency in the low tm region. It is a consequence of the larger Ht cut in Eq. ( 40 ) than 



in Eq. Q, which severely reduces the fraction of tagged events when mx ^ 300 GeV. 
Insofar as the Ht cut is present mainly to limit trigger bandwidth, one could also imagine 
lowering it and then controlling bandwidth by means of trigger-level 6-tagging. 
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Figure 11: Comparison bet-ween the hadron-level tagging efficiencies, at the 8 TeV LHC, for 
the toy Monte Carlo events and the radion and graviton MadGraphB events. We sho-w the total 
efficiencies and the breakup in different tagged samples, as a function of the boost factor tm- 



4.2 B tagging 

The final state that we are interested in includes four 6-quarks from the decays of the two 
Higgs bosons. Therefore, 6-tagging will be an important ingredient to improve the signal 
over background ratio. We have adopted in this study a 6-tagging scenario that we expect 



to be realistic (possibly conservative), inspired by the ATLAS and CMS capabilities 19 



20 87,88 . The probability of tagging a &-quark is taken to be fb = 0.75, the mistag 
probabilities of c-quarks, fc = 0.10, and of light quarks and gluons, fi = 0.03. We apply the 
6-tagging conditions on the parton level events after showering but before hadronization, 
that is, we tag b quarks rather than B hadrons. We will require one 6-tag per Higgs 
candidate. In detail it is implemented as follows: 

• Determine the number of 6-quarks within each of the two Higgs candidates' jets. 
Such candidate jets can be a single anti-/cT jet with radius R (in the boosted regime) 
or a jet composed by the sum of two different anti-A;^ jets (in the resolved limit). 
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Figure 12: The fo-tagging efRciencies for graviton mediated Higgs pair production at the LHC 8 
TeV. We show the efficiency of the resonance tagging (same as in the upper right plot of Fig. 11 1 
together with that including the 5-tagging. 



• A Higgs candidate jet is considered to be 5-tagged if it contains at least one b quark 
with pT^b ^ Px^b ~ GeV. The 6-tag efficiency is denoted by /;,. 

• A Higgs candidate jet which does not fulfill the previous condition, but contains at 
least one c quark with pT^c ^ Px^t^ ^^^^ 6-tagged with a mistag probability fe- 
rn A Higgs candidate jet which contains only light quarks and gluons will be 6-tagged 

with a mistag probability 

• 6-tagged events are those for which the two Higgs candidates' jets have been both 
6-tagged. 

Therefore, events will be given different weights according to the number of b and c quarks 
present in each of the two Higgs candidate jets. For instance, if the two Higgs candidates' 
jets each contain at least one b quark, the event is assigned a weight /^=0.56. This is the 
same probability for signal events and for QCD background events where two b quarks 
end up each in a Higgs candidate jet. 

We could also have considered a more optimistic scenario for the 6-tagging, in which 
each Higgs candidate is required to have two 6-tags. In particular CMS has demonstrated 
the ability to tag pairs of i3-hadrons even for angular separations AR^^i < 0.4 f87|, which 
suggests that this scenario could be viable also in the highly boosted limit where the 
S-hadrons are within a single anti-Zcr jet. 

Results for the 6-tagging efficiencies for graviton mediated Higgs pair production at the 
LHC 8 TeV as a function of rjv/ are shown in Fig. [12] As we can see, for the relevant mass 
range we have approximately a 15% constant signal efficiency after taking into account 
the resonance tagging algorithm and the 6-tagging. 

Strictly speaking, since there are two jets in a "resolved" Higgs candidate, the mistag probability 
is closer to 2/;. However, at the level of factors of two, our 6-tagging estimates for the backgrounds 
are probably not accurate. On one hand, for example, ATLAS [19[|20| obtains somewhat better light-jet 
rejection than the /; = 0.03 that we use. On the other hand, light-jet rejection factors will anyway depend 
on pt of the jets and potentially also their proximity to other jets. 
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QCD Dijets - LHC 14 TeV - PythiaS QCD multijets mistag probability - LHC 14 TeV 
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Figure 13: Left plot: the QCD dijet cross section before and after the events are processed 
through the resonance-tagging algorithm, for LHC 14 TeV. Right plot: the mistag probability of 
QCD dijet events with the resonance tagging (RT) algorithm, defined as the ratio of mistags over 
the QCD cross section, without and with 6-tagging, as a function of the mass point M . The mistag 
probability of QCD dijets at LHC 8 TeV is very similar. 



4.3 QCD multijet background simulation 

The dominant Standard Model background to multijet final states that leads to event 
topologies similar to the signal is QCD jet production. We have therefore produced a large 
sample of QCD multijets with PythiaS f44], starting from dijet configurations and with the 
shower radiation taking care of generating the higher-order jet topologies. The resulting 
hadron level events are then processed through the same analysis chain as the signal events. 
There are several ways in which QCD radiation can mimic the conditions for resonance 
tagging: for example, fake mass drops can be generated from a sufficiently symmetric 
splitting of a quark or gluon. Note that while PythiaS is known to underestimate the 
amount of QCD multi-jet topologies by a factor up to two with respect to experimental 
data [4], for the accuracy requested for this feasibility study we consider this precision 
to be sufficient. Similar results have been obtained with the Alpgen parton level event 
generator |89| matched to PythiaS using the MLM matching |90j. 



In Fig. 13 (left) we show the QCD dijet cross section obtained from the PythiaS multi- 



jet sample at LHC 14 TeV. The dijet cross section is defined, for each mass point M, as the 



number of QCD events that survive the basic selection cuts Eq. (40) and lead to an invari- 
ant mass within the mass resolution window around M given by [M{1 — fm), M(l + fm)]- 
In addition, we demand that the two leading jets are separated in rapidity by less than 
Aj/max- Note that the dijet cross sections flattens at small masses because there the 



selection cuts Eq. (40) have a sizable effect. 

In order to achieve an efficient QCD multijet event generation, for any candidate 
resonance mass M we have generated dijet events with a generation cut of pT > M/5, and 
no generation cut in rapidity. To motivate this choice, let us recall that the kinematics 
of massless jet pair production determine that the mass of the dijet will be given in 
terms of the pT of the jets and their rapidity separation in the laboratory frame Ay by 
M = cosh (Ay/2). Therefore, to properly cover all phase space the generation cut for 
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QCD dijets should be at least 



M 1 

2 cosh (Ay^ax/2) ^ ^ 

for any candidate resonance mass M. For the four jet configuration, it is reasonable to 
require a minimum pj- value of half of that of above. Since we are using a rapidity cut of 
^ymax = 1-3, we find that the PythiaS minimum pT in generation should be p™™ ~ M/5. 
We have explicitely verified that the QCD dijet cross-section is not modified if looser 
generation cuts are adopted. 



We also show in Fig. 13 (right) the background rejection factors, defined as the fraction 
of the QCD dijet events which are mistagged as arising from a heavy resonance, both with 
and without 6-tagging. Note that the background rejection probability is approximately 
scale invariant: similar mistag probabilities are obtained for all values of the mass. It is 
clear that the QCD background cross sections is reduced by a combination of the resonance 
tagging and 6-tagging by several orders of magnitude. For example, with 6-tagging the 
mistag probability is about 10~^, constant to very good approximation in all the relevant 
mass range. This improvement is due to the requirement that each Higgs candidate 
should be associated with two identified b-hadrons, a topology that is less frequent in 
QCD multijets. 

In summary, our study of the QCD background rejection confirm the consistency of 
the resonance tagging algorithm, since it makes possible to simultaneously explore the low 
mass and high mass region, achieving similar signal efficiencies and background rejection 
factors in all the mass range. 

4.4 Model independent exclusion limits 

We will now combine the results of the signal efficiencies and the multijet background 
rejection of the resonance tagging algorithm to derive model independent bounds on 
BSM scenarios with enhanced Higgs pair production in the 46 final state. This in- 
formation is enough to derive the values of the cross section times branching fraction 
a {pp — )• X) BR {X — )• HH) that can be excluded at the 95% confidence level from a mea- 
surement of the QCD 6-tagged multijet cross sections, as a function of the mass of this 
hypothetical resonance. In the following, to compute the number of signal and background 
events, we will assume a total integrated luminosity of £ = 25 fb~^ at 8 TeV and of 500 
fb-i at 14 TeV. 

For each candidate resonance mass, M, we compute the number of background events 
in a mass window of width fm =15% around M. The local p-value for each mass point M 
based on the expected number of signal and background events, Ns and N[, respectively, 
in the mass window considered, is given by 



1 

p 





[ 1 


^1 - Erf 





(42) 



2 

where Erf is the error function, and one assumes that the number of background events 



in each mass bin is a Poisson distribution with mean A*";,}^ Then requiring the condition 



Eq. ( 42 I is only valid where both Ns and Ni, are much larger than one, in the opposite case one has 
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p = 0.05 determines the number of signal events Ng that would allow an exclusion of the 
background-only hypothesis at the 95% confidence level, namely 



A^, = ^/Wk ■ Erf-i (1 - 2 • 0.05) . (44) 



Using Eq. (44) to determined the value of Ng in a given mass window, we can obtain the 
model independent bound on the combination a {pp — t- X) BR {X — )■ HH) by correcting 
the number of events for the signal tagging efficiency, with and without 6-tagging, the 
Higgs to hb branching fraction and the assumed total integrated luminosity £, 

I. ^ X, BR ,X ^ HH)U^,^ ,M) . iBmH^m^S^..mHM)C ' '''^ 

where the signal efficiency SignalEff (M) is derived from the MadGraph radion and graviton 



samples, see Fig. 12, We have used BR (if — )• 66) = 0.577 from the Higgs Cross Section 
Working Group recommendations ["ol]. 

The 95% excluded model-independent cross sections times branching fractions are 



shown in Fig. 14, We see that that we are sensitive to cross sections as small as 200 fb 
(50 fb) at M ~ 500 GeV at LHC 8 TeV (14 TeV), while at higher masses, M ~ 2 TeV, 
the 46 final state is sensitive to cross sections as small as 1 fb at both energies. Note 
that the increase in luminosity when going from 8 to 14 TeV is partially canceled by the 
corresponding increase of the high mass QCD dijet cross sections. On the other hand, 
signal cross sections in relevant models also increase when going from 8 to 14 TeV, so all 
in all we obtain a substantial improvement in exclusion reach when increasing the center 
of mass energy. 

These results confirm that the 46 final state can be relevant for many new physics 
scenarios that lead to enhanced cross sections for resonant Higgs pair production, and the 
search strategy that we propose makes it possible to efficiently explore a wide range of 
resonance masses within a common analysis. 

4.5 Graviton and radion searches in the 2H — t- 46 channel 

Now we consider the specific benchmark scenarios for radion and graviton production 
introduced in Sect. [Sj First of all, we summarize in Table [2] the model parameters that 
we adopt here. The mass scales and branching fractions to Higgs boson pairs are kept 
fixed, and only the couplings of the radion and graviton to gluons will be varied. For 
the graviton, we consider two different values of the coupling, Cg = 1 (G-Brane) which 
corresponds to the RSI model and Cg = 1/35 (G-Bulk) as in the bulk models. 

For the radion we study the nominal coupling Kg (R-Bulk) , as well as the case in which 
this coupling is enhanced by a factor of ten by some unspecified mechanism, such as when 
the radion arises as a composite bound state (R-Comp). We use the same mass scale in 
the two cases, = Ac = 2 TeV: although the two scales are related, we prefer to explore 



to use the correspong discrete Poisson formula for the p- value, 



which involves the incomplete Gamma function. 
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LHC 8 TeV, L = 25 fb ' 



LHC 14 TeV, L = 500 fb ' 




Candidate Resonance Mass ( TeV ) Candidate Resonance Mass ( TeV ) 



Figure 14: The model independent exclusion ranges in a {pp — )■ X) BR {X — ^ HH) at the 95% 
confidence level for the production of a heavy resonance X which then decays into a Higgs boson 
pair, which decay subsequently into a bottom-antibottom pair. The left plot show the exclusion 
ranges at 8 TeV with C — 25 fb~^ while the right plot corresponds to 14 TeV with C = 500 fb""*^. 
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Table 2: Parameters of the benchmark scenarios for radion and graviton production. For the 
radion we consider both the nominal value of (denoted by R-Bulk), and a coupling ten times 
larger that could arise for example in composite dual scenarios (denoted by R-Comp). For the 
graviton we consider two different values of the gluon-gluon-graviton coupling, Cg = 1 as in RSI 
(denoted by G-Brane settings) and Cg — 0.02 as in bulk models (denoted by G-Bulk). 



independently the radion and graviton scenarios. Let us recall that the cross section scale 
as so any different choice of the mass scale A.^,g will lead to a trivial rescaling of 

the cross section. 

First of all we evaluate the expected number of events for these four benchmark points 
at LHC 8 and 14 TeV, using the results of Sect. [3j We take into account the branching 
fraction of the Higgs bosons into bb pairs. As before, we assume total integrated luminosi- 
ties of £ = 25 fb"^ at 8 TeV and of 500 fb"^ at 14 TeV. The model cross sections for the 
benchmark scenarios can be easily obtained from the results of Sect. [3j in particular from 
Fig. [9} The number of expected events is shown in Fig. [TsJ after accounting for the selec- 
tion efficiencies from resonance and 6-tagging analysis. At 8 TeV we expect just a handful 
of events at low masses for the R-Bulk and G-Bulk points, and about one thousand events 
(few tens of events) and low (high) masses for the R-Comp and G-Brane points. At LHC 
14 TeV on the other hand we have a large enough number of events for all masses and all 
benchmark points, thanks to both the increased resonance production cross sections and 
the higher integrated luminosity. 

Now in Fig. 16 we show the same excluded cross sections at the 95% confidence level 
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LHC 8 TeV, L = 25 fb^ 



LHC 14 TeV, L = 500 fb ' 
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Figure 15: The expected number of events in tlie benchmark scenarios after both resonance tagging 
and 6-tagging, at the LHC 8 (left plot) and 14 TeV (right plot). The integrated luminosities are 
£ = 30 ft.-i at 8 TeV and L = 500 ft)-i at 14 TeV. 
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Figure 16: The 95% confidence level exclusion ranges at 8 TeV (left plot) and 14 TeV (right plot) 
compared to the specific cross sections of the four different model scenarios of Table [2] 



of Fig . [T4] but this time adding the specific model cross sections of the four scenarios of 
Table ^ The improvement in exclusion power when going from 8 to 14 TeV is clearp^ At 
LHC 8 TeV we can explore a large part of the parameter range of the graviton models 
up to 2 TeV, but the default radion scenario seems to be out of reach, unless its cross 
section is enhanced by some mechanism, for instance as in the composite duals discussed 
in Sect. IH At the LHC 14 TeV on the hand we are sensitive to R-Bulk scenario with 
Kg = 1, for most of the mass range up to Mfj) = 2 TeV. Likewise, we could exclude a Bulk 
graviton up to masses of 2.5 TeV. Therefore, after the energy increase to 14 TeV most of 
the parameter space of the radion and massive KK graviton models will become accessible 
in the 46 final state. 

Using these results, it is also possible to determine the 95% confidence level exclusion 
ranges for some of the parameters of the benchmark scenarios. We can keep all the 
parameters as in Table [2] and determine the exclusion ranges for the couplings of the 

^^Note that we only consider statistical errors in determining the exclusion limits. In the small and 
intermediate mass regions, the exclusions are based on a large number of events, corresponding to small 



signal over background ratios. This can be understood from Eq. ( 44 1 , which tells us that at the exclusion 
limit, Ns/Nb 



1/VNb 



30 



Graviton, 95% CL Exclusion Limits 



Radion, 95% CL Exclusion Limits 



10e 



10-' 





LHC 8 TeV 




RSI Cg=1 


LHC 14 TeV 


















Bulkc, 


=1/35 


3x10' 4x10' 


1 2 





Radion Composite |k,|=0.4 


LHC 8 TeV 




LHC 14 TeV 













3x10"' 4x10"' 



Graviton Mass ( TeV ) 



Radion Mass ( teV ) 



Figure 17: Left plot: the values of the massive KK graviton-gluon coupling Cg that can be excluded 
at the 95% CL as a function of the graviton mass for 8 and 14 TeV. For 14 TeV we provide the 
results for both fe-tagging scenarios. Right plot: same for the radion-gluon coupling k^. 



gluons to the massive KK graviton Cg and to the radion Kg and scan the allowed values 
for A^. 



We show the results in Fig. 17 In the case of the graviton coupling, we see that at 14 
TeV the 2H — )• 46 final state can access essentially all the relevant range, from the RSI 
value Cg = 1 down to the bulk value of = 1/35. For the case of the radion coupling Kg, 
we see that at 8 TeV we are sensitive to values down to Kg ~ 0.06 around 750 GeV, while 
at 14 TeV the LHC can exclude a bulk radion (with the default value for the coupling) 
for masses between 300 GeV and 2 TeV at least. 

In summary, we have shown in this section that the 2H — )■ 46 final state offers a 
promising channel to probe enhanced Higgs pair production at the LHC, despite the over- 
whelming QCD multijet background. The combination of jet substructure techniques and 
6-tagging makes it possible to probe a wide region in the parameter space of various bench- 
mark models. Therefore, we advocate that the experiments explore this new channel in 
order to complement existing searches of new heavy resonances in other, more traditional, 
channels. 



5 Conclusions and outlook 

In this paper we have presented a new strategy for heavy-resonance searches in multijet 
final states, which attempts to unify in a single approach the techniques used in the 
boosted and resolved regimes. By classifying events as a function of the number of mass- 
drop tags, we can smoothly interpolate between the boosted regime, where jet substructure 
techniques can be used, and the resolved regime, where the final state particles appear 
as well separated jets. In particular, we have considered the process X — )• YY — )• 4z in 
which the resonances Y are pair produced from the decay of a heavier resonance X and 
then decay into a pair of QCD partons, then observed as jets. We have shown that our 
strategy leads to approximately scale-invariant signal selection efficiencies and background 
rejection rates. 

As a benchmark scenario, we have considered Higgs pair production in extra dimension 
models, where the Higgs pair is produced from the decay of a heavy graviton or radion. 
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and then decays into four b quarks. Note however that the kinematical structure of the 
final state of the benchmark model holds for other scenarios, such as composite models, 
with more freedom on the couplings and therefore in the cross sections strengths. By 
comparing with the QCD multijet background, we have derived the model independent 
95% confidence level exclusion ranges for the cross sections for a {pp — t- X — )• HH) where 
the Higgs bosons decay into a bb pair, and showed that a substantial region of the parameter 
space of these models can be successfully explored in this final state with the tagging 
strategy that has been proposed. 

In the particular case of graviton and radion production, it would be especially interest- 
ing to study the feasibility of radion/graviton searches in the 6677 decay channel 
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This is a cleaner channel than the all hadronically decaying case, since the two high pT 
photons substantially reduce the QCD background. In this final state, by varying rM 
one also moves from the boosted regime (with one single fat jet in the final state) to the 
resolved limit, with two well separated jets in the final state. Another interesting final 
state to search for enhanced Higgs pair production would be bbZZ. 

The approach advocated in this paper could also be applied to other relevant problems, 
for example top quark pair production, again providing a smooth coverage acrosss the 
transition between the resolved regime, relevant for SM measurements, and the boosted 



regime, where substructure techniques [94 95 are used to enhance the potential of new 
physics searches. 
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